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Abstract 

We use an s-channel picture of hard hadronic collisions to investigate the parton distribution 
function for quarks at small momentum fraction x, which corresponds to very high energy scatter- 
ing. We study the renormalized quark distribution at one loop in this approach. In the high-energy 
picture, the quark distribution function is expressed in terms of a Wilson-line correlator that rep- 
resents the cross section for a color dipole to scatter from the proton. We model this Wilson-line 
correlator in a saturation model. We relate this representation of the quark distribution function to 
the corresponding representation of the structure function Ft(x, Q 2 ) for deeply inelastic scattering. 



1 



I. INTRODUCTION 



Proton structure at large Q 2 and small Bjorken x has been extensively investigated in 
experiments at HERA. This program is of great intrinsic interest and provides valuable 
information for the LHC program, where the short-distance structure of protons and nuclei 
will be probed at TeV energies. Two physical pictures that seem very different from each 
other are used to analyze hadronic structure functions for large Q 2 and small x. 

There is a parton picture, in which the hadron consists of partons and the partons undergo 
a hard collision that produces the final state. This is reviewed in [1]. This applies at large 
Q 2 . The corresponding theoretical method is that of factorization. The cross section is 
written as a convolution of parton distribution functions, f a /A(x,Q 2 ), and a hard cross 
section for partonic scattering. The parton distribution functions are evaluated at small x, 
but the x dependence is not predicted except insofar as it results from evolution starting 
from fa/A^x, Ql) at a smaller virtuality scale Q^. 

There is an s-channel picture, in which one thinks of the event in the rest frame of one of 
the hadrons. See [2] for recent accounts. This applies at small x. The hard interaction takes 
place far outside the hadron and the products of the interaction travel toward the hadron 
and interact with it. In the simplest case, there are effectively two objects that collide with 
the hadron. These objects carry opposite color, so that they can be said to constitute a 
"color dipole" , described by a correlator of two eikonal Wilson lines. An important concept 
here is that the cross section for the color dipole to scatter from the hadron can be simple 
when the transverse separation between the elements of the dipole are large. Then the dipole 
always scatters as long as its impact parameter is within the hadron radius. One speaks of 
the cross section saturating - that is, being as large as it possibly could be [3, 4]. 

These pictures seem quite different, as they look at the collision in different reference 
frames, and lead to different theoretical methods, but they are not at all incompatible. In 
the region where their domains of validity overlap, they must describe the same physics. 
The aim of this paper is to connect the two pictures. 

We examine one of the main ingredients used in the parton picture, namely the distri- 
bution function for finding a quark in a hadron, defined as a hadronic matrix element of 
a certain product of operators [5]. For very small x, the parton system created by this 
operator is far outside the hadron. We analyze the evolution of the system using s-channel 
methods. We find that the quark distribution can be expressed as a Wilson-line correlator 
convoluted with a simple lightcone wave function. Moreover, we find that this answer allows 
one to relate with precision the seemingly dissimilar results for structure functions in the 
parton framework and the s-channel framework. 

Part of the results of this analysis have been used in [6] to investigate the power corrections 
to structure functions that arise from the s-channel picture. 

The content of the paper is as follows. We begin by applying the hamiltonian method [7] 
to the quark distribution function. This allows us to write the quark distribution as a 
convolution of a lightcone wave function and a matrix element of eikonal-line operators 
(Sec. II). We work in the lowest-order approximation, i.e., the dipole approximation. The 
convolution formula provides a simple interpretation in coordinate space for the physical 
process that probes the distribution. The parton distribution is defined by matrix elements 
of operator products that require renormalization. We perform the renormalization at one 
loop using the MS subtraction scheme for the ultraviolet divergences. 

The eikonal-operator matrix element receives contribution from both short distances and 
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long distances. We first analyze it by an expansion in powers of gA (with g the strong 
coupling and A the gauge field), valid at short distances (Sec. III). This expansion is useful 
to carry out the matching with renormalization-group evolution equations. In particular, it 
allows us to relate the eikonal matrix element at short distances to a well-prescribed integral 
of the gluon distribution function. 

Next, we motivate and discuss a widely used approach for modeling the eikonal-operator 
matrix element at large distances (Sec. IV), based on parton saturation [3, 4]. As the 
saturation scale in the quark sector is likely to be at much lower momenta than in the 
gluon sector (see e.g. [8]), we critically examine the validity of the treatment for the quark 
distribution, and the potential breakdown of the dipole approximation (Sec. V). 

We finally discuss the relation of our results for the quark distribution with known dipole 
results for structure functions (Sec. VI). This discussion also illustrates how standard fac- 
torization properties are reobtained from the s-channel point of view. 

Some supplementary material is left to the appendices. In Appendix A we collect calcula- 
tional details on integrals of lightcone wave functions. In Appendix B we give a relation be- 
tween products of eikonal operators for color-octet and color-triplet dipoles. In Appendix C 
we report details on applying the hamiltonian method to the hadronic matrix element of 
two currents. 



II. QUARK DISTRIBUTION IN THE S-CHANNEL PICTURE 

We study the quark distribution using the s-channel picture in the style of [7]. We start 
with the definition [5] of the quark distribution as a proton matrix element of a certain 
operator, 

x (P, S |^(0)Q(0) T +Q %~)V<0, y~, 0) \P, s) c . (1) 
Here for any four-vector z M we use lightcone components z^ defined as 

The proton momentum is 

/ M 2 \ 

P=(P + ,P-,Pj = [P + ,^0) ■ (3) 

The operator is the path-ordered exponential of the color potential 

QV)=pexp{-^ J^°° dz-A+(0,z-,0)t o } , (4) 

where the path ordering instruction V puts fields and color matrices with the most positive 
values of z~ to the left. Equivalently, following the notation of [7], we can think of Q^{y~) as 
creating an eikonal particle that moves in the minus direction, starting at minus coordinate 
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y~ . An eikonal particle is an imaginary particle that retains its plus and transverse positions 
no matter how much momentum it absorbs. The subscript c on the matrix element in Eq. (1) 
indicates that we are to take the connected parts of the graphs, in which some partons from 
the proton states communicate with the indicated operators. The operator product in 
Eq. (1) is ultraviolet divergent and requires renormalization. We will use the standard MS 
prescription. The required subtraction at the one loop level is analyzed in Sec. HE. 



A. The quark distribution as a forward scattering amplitude 

We begin by rewriting the matrix element in Eq. (1) so that it has the form of the real 
part of a forward scattering amplitude. To do this, we write f q / p (x, /S) in two pieces, 



where 



and 



We note that 
Thus / is twice the real part of / 



f q/p (x, fj) = f+ /p (x, n) + f q/p (x, fj), (5) 

x (P, s|<K0)Q(0) 7 + Q VMo, y~, 0)\P, s) c (6) 

x (p,s\mQ(oh + QKyim v~, o) \p, s) c . (7) 

fq/p( X ^) = ifg/ p ( X ^)T- (8) 



x(P, S |T{^(0)g(0) 7 + g t (r)^(0,r,0)} \P,s) c . (9) 

The operator product in /~ is time ordered since y~ < 0. The T here indicates this time 
ordering. For our purposes, it is helpful to insert a factor x and another y~ integral: 



(P', s\T{^(0, y;, O)g( 2 / 2 -) 7 + g t (2/r)^(0, yT, 0) } \P, s) c . (10) 



Re 



The integral over P' + can be thought of as setting the proton state to position y~ = 0. 
We have thus rewritten the original /, which was analogous to a total cross section, as a 
Green function analogous to a forward scattering amplitude. Our next task is to break the 
scattering amplitude into parts that can be analyzed separately. 
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B. Decomposition of the gluon field 



The Fourier transformed operator Q^(i/i)ijj(0, y^, 0) in Eq. (10) creates an antiquark and 
an eikonal particle with a total plus-momentum xP + . We consider that x is very small, say 
10~ 3 . That means that the typical distance \f[ from the proton to where the antiquark and 
the eikonal particle are created is large, of order l/(xP + ). This is way outside the proton. 
The antiquark and eikonal particle develop into a shower of partons with minus-momenta 
of order k~ = (k\ + k 2 )/(2k + ) ~ m 2 /(xP + ), where k± is the transverse momentum of the 
parton and k 2 is its virtuality and we take both of these to be of order m 2 = (300 MeV) 2 . 
Thus the partons created by the original operator have very large minus-momenta. We will 
speak of them as "fast" partons. As noted, the fast partons travel a long distance in y~ 
before meeting the proton. 

When the fast partons meet the proton, they scatter from the gluon field of the proton, as 
depicted in Fig. 1. The gluon field of the proton consists of "slow" gluons, with plus- momenta 
much larger than xP + . (Then the minus-momenta of these gluons, k~ = {k\ + k 2 )/(2k + ) is 
much smaller than m 2 / (xP + ), assuming again that k\ and k 2 are of order m 2 .) We represent 
the gluon field produced by the proton as an external field A^(x) and consider the quantity 

U[A] = ^ fdy^dy^e^ p+ ^^ 

x { (o|^(o, y^, o)g(y 2 -) 7 + g f (yl)m yT, o)\o) A 
- (o\m y 2 - oM^-h+Q^rMo, yr, o)|o> } . (n) 

This is the amplitude for the fast partons to be created by the operator Q'ij), scatter from 
the external field A, then be annihilated by the conjugate operator ipQ. In the second term, 
we subtract a no-scattering term with the external field set to zero, in accordance with the 
instruction to take only connected graphs. Then the quark distribution is a proton matrix 
element of U[A], with the external field A replaced by the quantum field A, 

There is, of course, a catch in this. There is only one gluon field A^(x). We need to divide 
it into two pieces, one associated with the fast partons and one associated with the proton. 
To do this, we choose a momentum fraction x c . Gluons with plus-momenta smaller than 
x c P + are associated with the fast partons. Gluons with plus-momenta larger than x c P + are 
associated with the proton and included in the external field A in Eq. (11). In order for 
the approximations discussed below to work, we need x x c . It is perhaps easiest to think 
about the physics taking x c to be of order 1/(R P P + ), so that the proton's field is considered 
to have a spatial extent of the order of the proton radius, R p . However, in the end we will 
want to take i c «l and in fact let x c be pretty close to x. 

We will, in fact, not need to be very specific about how to implement the division at 
momentum fraction x c . For our discussion of the quark distribution, working at lowest 
order in perturbation theory, we are saved from sensitivity to the splitting method by the 
fact that the g — > qq Altarelli-Parisi splitting function does not have a soft singularity. If we 
worked with the gluon distribution or with the quark distribution to higher order, we would 
need a more sophisticated analysis. 
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FIG. 1: Quark distribution in the s-channel picture. The black circles at the top represent the 
insertion of the operators in Eq. (10). The double straight line represents the eikonal. Any number 
of slow gluons couple the fast parton and the eikonal to the proton. 

C. The evolution operator U at high energy 

The function U[A] can be written in the interaction picture with A as the perturbation: 

u[A] = — J dtadytefa - yi) e-^for-in 

x{0\U(oo,y^(0,y^,0)Q(y^+ 

x E% 2 - yr)Q + (yrM0, Vl, 0)U(yI, oo)|0) 

— (same with A= 0). (13) 

This is without approximation. Now we recognize that for small x, only y^ ^> r~ and 
y^ <C —r~ are important, where r~ is the effective radius of the proton's field in the 
longitudinal direction, r~ = l/(x c P + ). The external field is concentrated in \y~\ < r~, so 
we have 

U[A] « ^- / dy~ 2 / dyle-t^-^ 

^ JO J-oo 

x (0^(0, y 2 -, 0)g(y 2 -) 7 + [f/(oo, -oo) - l]Q^y^(0, y^, 0)|0). (14) 

Here we have subtracted the no-scattering term as "—1". 

At this stage, our quantum fields are evolving with full QCD [not including the external 
field, which is represented in U(oo, — oo)]. Let us now expand this evolution in powers of a s 
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and take just the Born term. Then the fields evolve with just the free field hamiltonian. We 
can insert intermediate states, and at this level of approximation, the intermediate states 
contain just one antiquark and the one eikonal particle £. We get 

U[A] » ^- J d y- 2 j dy^e-™ p+ ^-y^ 

x <O|-0(O, y 2 , 0)Q(y 2 h + \p 2 , p 2 , s 2 , 8} 

x{p 2 ,P2,s 2 ,S\U(oo,-oo) - l|p7,pi,si,£) 

x (pi , Pl , Sl ,E\Qi(y^(0, yr, 0) |0). (15) 

Taking into account that particle 1 has plus momentum pf = p\/[2p^) while particle 2 has 
plus momentum p 2 = p 2 / (2p 2 ), we can evaluate the dependence of the matrix elements on 
and y 2 as 

/*oo /*0 

x / dy 2 / ^ re - J ( a;P++ ^)^e +j(a;P++ ^ ) ^ 

JO J-oo 

x(o|^(o)g(o) 7 + b2,P2,s 2 ,^) 

x (P 2 ',f»2,S2,^|f/(00,-00) - l|p]~,Pl,Si,£) 

x^Pi^i^IQ^OMO^O). (16) 
We can now perform the y~ integrations to produce energy denominators: 



—i —i 



x 



xP+ + p 2 xP + + pj~ 

x(o|^(o)g(o) 7 + |^,p 2 ,s 2 ,^) 

X(P2 >P2, S 2 ,^|f/(00, -OO) - l|pf,Pi, Si, 5) 

x{pr,pi,5i,5|Qt(0)V(0)|0). (17) 
For the factor giving the interaction of the partons with the external field, we have 

(P2">P2,S2,£|^(00,-00) - llp^PljS!,^) = 



where 



and 



2n2p^5(p^ -p 2 )5 SlS2 [F( Pl -p 2 ) t F(0) - (2tt) 2 5(Pi -p 2 )], (18) 
F(A) =Vexp^-igJ + dz~A^(0, z~, A)t J (19) 
F(fe) = / (IAe* A F(A). (20) 



This gives 

U[A] « J dA -L Tr[l - F(A)tF(0)] u(A), (21) 



where 



M(A) = w jT dp ~ / rfP2 / dpi £ e * A ' (pl ~ P2) 



Thus 



where 



(2xP+p~ + p2)(2a;P+p- + pf) 
x(0|^(0)g(0) 7 + b-,p 2 , s,£)(p- pi,s,£|Qt(0)V>(0)|0). (22) 

xf q/p {x,n) = ReJdA u{A) 5j(A), (23) 



s ^)= (sEj / (^^4 Tr[l - F(A)tF(o)]|p ' s) - (24) 

Here F(A) is now defined with the quantum field A, 

F(A)=Vexp{-ig J + dz~A+(0, z~, A)t tt J . (25) 

Eq. (23) has a simple interpretation. First, m(A) is the square of the antiquark wave 
function, giving the probability that the antiquark has reached a separation A from the 
eikonal line by the time it reaches the hadron. Second, we have a probability Sj(A) for the 
antiquark-eikonal dipole to scatter from the proton. 

D. The squared wave function for the antiquark 

We now need the function u(A). First, we need the operator matrix elements: 

m = ^(o|^(o)Q(o) 7 + b-,p 2 , s ,f)(p-,i> 1 , s ,£:|gt(o)^(o)|o) 

s 

= ^(0|^(0)7 + |p-,P2,s)(p-, Pl , S |^(0)|0). (26) 

s 

There is an implicit color trace here. Restoring the color indices makes it 

M = J2( Q \M0)l + \p~,P2,s,P)(j ) -,p 1 ,s,P\M0)\0). (27) 

s 

Writing this with spinors gives 

M = S a(i 5p a ^2v(p~ ,p 2 , s)j + v(p~ ,pi, s). (28) 



8 



Now we need to know about the spin states. We use null-plane helicity states appropriate 
to the the "time." These have the normalization 



Thus 



v(p ,£»2,s')7 v(p ,p 1 ,8)=2p 5 SS >. (29) 



M = 5 a p6j3 a ^2v(p ,p 2 ,s)^ + v(p ,pi,s')5 ss > 

ss' 

= 7T-/v(P ,P2,s)-y + v(p ,p 1 ,s')v(p ,p 2 ,s')~/ v(p ,pi,s) 

ss' 

6 Tr{^7 j) 2 l + } 



2p 



2p~ 

Pi ■ P2- (30) 

p- 



Thus 



u(A) = -j—^- J dp~jdp 2 jdp 1 e i ^-^ 



P 



(2xP+p~ + pj)(2xP + p~ + p\) 

2N C 
x — p 1 ■ p 2 

V 

AN roc 



POO p p 

I dA 2 J dp 2 J d Pl e iA < P1 -^ 



{x + pl){x + piy 



where we have defined A 2 = 2xP + p . Extending this to 4 — 2e dimensions, we have 



_ an c ^ r°° r r A . (P1 _ P2) 

U(A) " (2k)*-± h J P2 J Pl (A 2 + Pi) (A 2 + Pi) 

We perform the integration separately in Appendix A. We find 

^ c i / ,. 2A2 , 2e r(2- e ) 2 



u(A) 



( T ^AT 1 7 . (33) 
37T 4 A 4 V ' 1 - 2e/3 V ; 



E. Renormalization of the quark distribution 

Using the result (33) for u(A), we have 



= 3^ d 2 -^ i (vr/, 2 A 2 ) 2e S 7 (A) - UV. (34) 
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Here there is an ultraviolet divergence from the small A integration region. The notation 
indicates that we should renormalize the divergence by subtracting a UV counter term. 

The standard definition of the parton distribution functions gives these functions as 
hadron matrix elements of operator products that must be renormalized [5]. It is thus not a 
surprise that we have a divergent integral in Eq. (34). The standard treatment is to apply 
MS renormalization. At the one loop level at which we work here, this means performing 
the integrals in 4 — 2e dimensions and subtracting a counter-term of the form 

"UV" = const, x - ^ >_ . (35) 

In this section, we implement this subtraction, turning it (approximately) into a cutoff on 
|A|. 

We are eliminating only the divergence from the innermost loop in the Feynman diagrams 
that define xf q / p (x, /i), so we treat the outer loops in Sj(A) as containing only soft momenta. 
For this reason, we treat H/(A) as being an analytic function of A near A = 0. We thus 
write 



lj(A) = S,(0) + A t [^S / (A)] A=0 + iAW[^^H J (A)] A= o 

2 (36) 



+ lAWA fc [^ fc H 7 (A)] A=0 + R(A). 

(We follow the convention that indices k are summed from 1 to 2 or, with dimensional 
regularization, from 1 to 2 — 2e.) The remainder, R(A), goes to zero like A 4 as A — > 0. The 
first term vanishes because 2/(0) = by construction. The second and fourth terms vanish 
upon integrating over A. In the third term, under the integration over A, we can replace 

A i A j -> — !— S' l ' J A 2 . (37) 
2 - 2e v ' 

Thus, in the small A integration region, we can replace 

S/(A) - j^—j A 2 [d 2 ± ~j(A)} A=0 + R(A). (38) 

We introduce this approximation in the small A integration region, defined by A/i < a, 
where a is a parameter of order 1 that we can adjust. Thus we write 

Nc T, 2a „ /a2 , Sj(A) 



N f 

xf q/p {x,v) = ^ d 2 A 9(A 2 fi 2 > a 

+ Jd 2 A9(A 2 fi 2 <a 



A 



i 



R(A) 



+ T ^[^(A)] A .o A ^ 



r(2 - e) 2 2e f 2 _ 2 6(A 2 ^ 2 < a 2 ) 

^ / d A A2 V*^ A ) 



(l-e)(l-2e/3) ^ J A 2 
- UV. 
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In the first two terms, there is no ultraviolet divergence, so we set e — > 0. We can perform 
the integration in the third term to obtain 



■ J d 2 A J(AV < «^ 



3vr 4 7 A 4 

r(2 - e ) 2 



e)(l-2e/3) 



(«74) £ - 1 



Here we have identified the UV subtraction term and written it as —1 inside the braces in 
the last term. The —1 removes the 1/e pole, but, in general, leaves a remainder that is finite 
as e — > 0. We set 

a = 2e 1/6 " 7 « 1.32657, (41) 

where 7 is the Euler constant, 7 ~ 0.577216, that appears in T(l — e) = 1 + 67 + • ■ ■ . With 
this choice, we cancel the finite term and leave 



xf q/p (x,») = ^ jd 2 A fl(AV > a 2 ) ^ 



(42) 



3tt 4 / v n A 4 



The term containing R(A) is needed to express the result of MS renormalization if is 
of the order of the proton radius, R p . However, the parton distribution function evaluated 
at such a renormalization scale is not really a very interesting object. For large values of 
/i, the term containing R(A) is of order l/(/i 2 i? 2 ) and can be neglected. Thus, as long as 
l/(n 2 Rp) <C 1, we can write 1 

xf q/p ( X , ,x) = ^ JdA £(AV > a 2 ) (43) 

This is a remarkably simple formula. Almost everything is contained in the dipole scattering 
function Sf(A). In the following section, we will study S/(A) for small A, where pertur- 
bation theory can be used. Then in Sec. IV, we will introduce and motivate on physical 
grounds a well known model for H/(A) for large A. 



III. DIPOLE SCATTERING AND THE GLUON DISTRIBUTION 

In this section, we investigate the dipole scattering function S/(A) for small A, where 
the use of perturbation theory is allowed. We will see that Hj(A) for small A is related to 
the gluon distribution. 



The derivation assumes that S/(A), and thus -R(A), is defined with a fixed renormalization scale fi. If 
we use Eq. (43), however, we can set \i in S/(A) to a/ A. 
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A. The dipole scattering function at small A 



We begin by studying H/(A) for small A and at lowest order in an expansion in powers 
of the strong coupling g. For the sake of generality, we consider using color matrices t c 
in a representation r of SU(N) that need not be the fundamental representation that is 
appropriate for the quark distribution. We start by defining SJ(A) corresponding to the 
representation r, 

(2tx)2P i+ 5(P + - P I+ )5 SS , cf r Hj(A) = (P\ S '|Tr[l - F r (A)tF r (0)]|P, s), (44) 

where + ^ 

^r(A) = Pexp <y—ig J dz~A+(0, z~, A)t a J . (45) 

The matrices t a here are in the color representation r and d r is the dimension of the repre- 
sentation, 

fundamental : d r = N c , (46) 

adjoint : d r = iV 2 - 1 . (47) 

For the fundamental representation, averaging over spins and integrating over P' gives the 
definition (24) of Sj(A) that we used for the quark distribution. 

We are interested in small A, for which an expansion in powers of gA is justified. If we 
limit ourselves to evaluating the trace to the accuracy {gA) 2 , we can ignore the P-products 
and the non-commutativity of the fields in the exponent: 



Tr[l -P r (A) f P r (0)] « Tr 



1 - exp \ ig j dz~ [A+(0, z~, A) - A+(0, z~,0)] t)j 



■ (48) 



The exponent has the expansion around A = 

A c + (0,z-,A)-A c + (0,z-,O) = A i [5 i A+(0^-,A)] A=0 

+ iA'A^ftfl^CO, 2", A)] A=0 + ■ ■ • . (49) 

Now we evaluate the trace to order A 2 by expanding the exponential. The zeroth order 
term from the exponential expansion cancels against 1. The linear term gives zero because 
Tr t c = 0. The quadratic term in the exponential expansion receives an order- A 2 contribution 
only from the first-derivative term in the exponent (49), again because Tr t c = 0. Thus 



Tr[l-F r (A)tF r (0)] ~ --(i 5 ) 2 A ! A^Tr[fA] / dzf / dz^ 

^ J— oo J —oo 

x [^+(0,zr,A)] A=0 [^A + (0,^,A)] A=0 . (50) 

The trace of two generators is 

Tr[t a t b ) = c r 5 ab , (51) 
where c r depends on the representation: 

fundamental : cv = Tr — 1/2, (52) 

adjoint : c r = Ca = N c . (53) 
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For matrix elements over hadron states carrying no transverse momentum we can replace 

AWd'A^At - l -A%A + a d,At. (54) 
Thus, to second order in gA and second order in A, we have 

(P', s'\Tr [1 - F r t(A)F r (0)] \P, s) « -g 2 c r A 2 / cfef / dzj (55) 

^ J —oo J— oo 

x (p', s '|a,-A+(o, *r, 0)^A+(0, 2T 2 -, 0)|P, S ). 

At this point we make explicit our restriction on the field operators A M (x), namely that 
only modes with gluon momenta larger than x c P + are included. (See Sec. II B.) This 
means that the reach in coordinate space in the integrations over and z% is limited to 
\zi — z^ \ < l/(x c P + ). Thus we write 

2 poo poo 

(2n)2P' + 5(P + - P' + )5 ss/ Hj(A) » A 2 / dz~ x \ dz~ 2 

4 Ci r J -oo J —oo 

x < V^cP" 1 ")) (56) 

x (P', s'l^A^O, *r, 0) ^A+(0, 0)|P, s). 



B. The gluon distribution function 

One might suspect that the right hand side of Eq. (56), being quadratic in the gluon 
field, may have something to do with the gluon distribution function. Indeed, as we shall 
see later, the relation is well known. We can check this relation by referring directly to the 
definition of the gluon distribution [5], 

fg/p&n) = 2~p+ J dy-e^- ( ± £j (P, s\F^(0, 0, 0)P>(0, y~, 0)|P, s), (57) 
with 

F a (y) +J = E{y) ab F b {y) + \ (58) 



E(y)=Vexp(-ig f dz~ Aj(y + , z~ ,y)t)j . 



(59) 



We study this distribution in the s-channel picture. At the lowest order in this picture the 
gluons from the background field couple to the vertex measured by the operator (57). 

We use momentum conservation to insert a second integral over the minus coordinate in 
Eq. (57), and we also use rotational invariance to eliminate the spin average: 

1 POO POO 

2n5(P' + - P + )5 ss ,f g/p (x, /i) = J dy^ J dy^ e * xP+ ^^ 

x (P', S '|P^(0,^,0)P^(0,yr,0)|P, S ). (60) 
As in [7], it is convenient to rewrite F +J as 

KKv) = d + (E(y) ab Al(y)) - E(y) ab d> A+(y) (61) 
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and note that inside the integral in Eq. (57), d + gives —ixP + . Thus, in the limit x — > 0, 
the first term in Eq. (61) can be neglected. Additionally, to lowest order in a perturbative 
expansion, the eikonal operator E(y) a b is equivalent to the unit operator. Thus we replace 

F+*(y) - -&At{y). (62) 

This gives 

2tt5(P /+ - P + )5 ss ,f g/p (x, /i) « J dy^ J dy 2 e ixP+ ^~^ 

x (P',s'\d j A+(y)(0,y^,O)d j At(y)(0,y^,O)\P,s). (63) 

We need one more approximation. For small x, the factor exp(ixP + (yi — y^)) is ap- 
proximately 1. This is not exact, and fails for very large (yjf — y^)- For \y^ — y% \ > 
l/(xP + ), the matrix element is a slowly varying function of (yf — yj)> so the oscil- 
lating factor exp(ixP + (y± — y^)) effectively cuts off the integral. Thus we approximate 
exp(ixP + (yi — y^)) by a theta function that restricts the integration to \y{ — y 2 \ < 
l/(xP + ). This gives 

/*oo poo 

27t5(p /+ -p + )8 ss ,f g/p (x,fi) « — — y dyry ^ 2 -^(br-y 2 -|<i/(xP + )) 

x (P\ s'\d,At(y)(0, y^ 0) a,v4+(y)(0, yf, 0)|P, s). (64) 



C. Relation between f g / p and E 

If we compare Eqs. (56) and (64), we see that 

ng 2 c r 2 



S /(A) ~ -^--j-A 2 x c f g/p (x cl Li). (65) 



Eq. (67) in particular implies that 

Hj,fund. ^ (l/2)/3 _ 4 CV 
5/,adj. 3/8 9 Ca 



(66) 



for small A and small x and to lowest perturbative order. There is a more general relation 
between H^nd. and S/ adj., which we give in Appendix B. 

The case of interest to us here is that of the fundamental representation, which applies 
to the quark distribution function. For this case, the result is 

- ^^x c f g/p (x c ,fi). (67) 



D. Matching using the renormalization group 

In deriving Eq. (67), we have employed rather crude approximations relating to the 
integrations over the minus component of position for the gluon field. The main idea was 
that the structure in matrix elements of A + {y) occurs for y~ less than l/(x c P + ), so that 
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limits on the integrations over y~ should not much matter. If this were precisely the case, 
then the function x c f g / p (x c , fi) that appears on the right hand side of Eq. (67) would be 
independent of x c . In fact, x c f g / p (x c , /i) grows slowly as x c decreases. Thus, we should try 
to make the relation (67) more precise. 

We can use the scale dependence of f g / p (x,fi) to provide a more precise matching condi- 
tion. On one hand, we have the leading order renormalization group equation, 

d ot f x fx \ a f x fx \ 

J^fMx^) = ^ I dzP ^- z f^{- z ^) + ^J x dzP qq (z)-f q/P {-^)- (68) 

At small x, the gluon distribution dominates and the quark distribution is effectively a s times 
the gluon distribution (as is, in fact, consistent with this equation). Thus the renormalization 
group equation can be approximated by 

where we have inserted the specific form of P qg {x/y). 

In our small x approximations, f q / p (x,fi) is given by Eq. (43) (as long as « ^> \/R p ). 
Differentiating this equation with respect to « gives 

Comparing these equations gives 

I^Ug^fea/A), (71) 

for A -C R p , where 

xG(x, = 1 J'dz [z 2 + (l- zf] % - f 9/p //) . (72) 

Note that the lower limit on the z integral is just a reminder that yf g / p (y,fi) vanishes for 
y > 1. Note also that the integral of the weight function is 

\f Q dz [z* + (l-z) 2 ]=l. (73) 

Thus xG(x,fi) is yfg/p(y,fi) averaged over values of y that are somewhat larger than x. If 
we consider a typical value of z to be 1/2, then the typical value of y at which the gluon 
distribution is evaluated is y = 2x. If, for example, xf g / p (x,/j,) oc x~°' 3 for small x, then 
xG(x,fi) ~ 0.76 x xf g / p (x,fi) for small x. 

Eq. (71) is the same result as in Eq. (67), except that now x c f g / p (x c , /x) is replaced by the 
more precise value, xG(x, a/A). Note that the matching condition suggests that x c be set 
to a value not much bigger than x. This is in part because our perturbative calculation of 
U (A) was to zeroth order only. Had we worked to one more order in perturbation theory, we 
could have included the emission of a fast gluon with momentum fraction between x and x c . 
However, to the order to which we calculated, there were no interactions with fast gluons. 
Working to this order, the best choice is to include all possible gluons as slow gluons. This 
means setting x c to something close to x. 



SKA) 
A 2 



A=a/ it. 



(70) 
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IV. THE HADRONIC MATRIX ELEMENT 

In this section, we motivate a widely used model for S/(A) that applies at large A. 
We begin by writing 5/ (A) as an integral of a function H(6, A) that has a direct physical 
interpretation. 



A. Scattering at fixed impact parameter 

We can rewrite the hadron matrix element in Eq. (24) by introducing an integration over 
an impact parameter 6. Denoting an eigenstate of transverse position by a subscript x, we 
have 

(P' + , 0, s\±- Tr[l - P(A)tp(0)]|P+, 0, s) 
= (P' + , 0, s\±r Tr[l - F(A/2)tF(-A/2)]|P + , 0, s) 
= ±- Jdb x (P f+ ,-b,s\Tr[l- F(A/2)^F(-A/2)]\P + ,0,s) 
= ±r Jdb x (P' + ,0,s\Tr[l-F(b + A/2yF(b-A/2)]\P + ,0,s). (74) 

dbZ(b,A), (75) 



Thus 



where 2 



S(6,A) = i-^j |^^^0, S |Tr[l-F(6 + A/2)tF(6-A/2)]|P+,0, S ). 

(76) 

The quantity 5(6, A) is more suitable than S/(A) as a quantity to model since the physics 
of the dipole-proton interaction should depend on b. Given a model for 5(6, A), one obtains 
5/ (A) by integrating over b. 

Given that 5/(A) has the behavior given by Eq. (71) at small A, we can write for 5(6, A) 
at small A, 

5(6, A) = A 2 ^ xG(x, a/A) 0(6), (77) 

where 

^6 0(6) = !. (78) 



Given that xG(x,a/A) is the number of gluons per unit d\ogx (averaged over momentum 
fractions somewhat larger than x), we interpret xG(x, a/ A) 0(6) as the number of gluons per 
unit area db and per unit d log x at a distance 6 from the center of the proton. Consistently 
with this interpretation, we assume that 

0(6) > (79) 



Here S(b, A) is the spin average (with s 1 = s) of what is called c in [8]. 
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and 



We will need a model for 0(b). 



0(b) = for |b| > R p 



(80) 



B. Interpretation and properties of 3(6, A) 
Let us write S(b, A) as 

S(6,A) = l-T(b,A). (81) 

Here the 1 comes from the 1 in Eq. (76). Then T comes from the matrix element of F^F. In 
the language of classical optics, T(b, A) is the transmission coefficient for a dipole of size A 
impinging on the proton at impact parameter b. According to the definition (76), the dipole 
is counted as transmitted only if the proton is left intact after the dipole moves through it. 
(This is the consequence of our having switched from a description of f q / p (x,fi) as a total 
cross section to a description in the form of a forward scattering amplitude.) Based on this 
interpretation and on what we have already learned about S(b, A), we expect S(b, A) to 
have the following properties. 

1. T(b, A) = 1 for \b\ > R p + A/2. 

2. T(b, A) = 1 for A = 0. 

3. T(b, A) m for \b\ < R p with \b\ not close to R p and A not small. 

4. T(b, A) = 1 - A 2 \K 2 a s /(2N c )} xG(x, a/A) 0(b) + C(A 4 ) for A -> 0. 

Property 1 simply says that a dipole that entirely misses the proton does not interact with 
it and is thus perfectly transmitted. Property 2 holds because a dipole with zero separation 
does not have any interaction with the proton. This is the property of color transparency. 
Property 3 applies because a big dipole has strong interactions, so that we expect that after 
such a dipole moves through the proton the proton is almost never left intact. Property 4 is 
consistent with T being 1 for A = and reflects our previously obtained perturbative result 
for H at small A. 



C. Model for 3(6, A) 

There is a simple model for T(6, A) that is consistent with the properties listed in the 
previous subsection, 

T(b, A) = exp (- A 2 ^ xG(x, a/A) 0(b)) . (82) 

This is a small variation on the widely used saturation model [4, 9], with the gluon dis- 
tribution treated according to the matching of Sec. HID. The same model for H(b, A) 
is 

2(b,A) = l-e- A2Q ^/ 4 , (83) 
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where Q s (b), known as the saturation scale, is 

Ql(b) = 2 -^xG(x,a/A)<J ) (b)- (84) 

This is the saturation scale for a dipole in the fundamental representation. From Eq. (66), 
we have for a dipole in the adjoint representation (as would be appropriate for the gluon 
distribution), 

Q 2 s (b, adjoint) = ^ Q 2 s (b, fundamental) = ^ xG(x, a/A) 0(6). (85) 

Of N c — 1 

The name of the model and of the scale Q s derives from the fact that S(6, A) grows as A 
increases until it saturates with S(b, A) pa 1 when A reaches approximately 2/Q s . 
For a specific model, we follow Mueller [4] in choosing 

m = 2^ " b2 m < Rp) - (86) 



V. CRITIQUE OF THE MODEL 

The dipole picture and saturation model [4, 9] along the lines just described has enjoyed 
some success when its predictions are compared to experimental results in both inclusive 
and diffractive deeply inelastic scattering ([3], and references therein). We do not attempt 
a numerical comparison in this paper. However, we do offer some comments on the extent 
to which the dipole picture for f q / p should be expected to be reliable. 

We have found that the parton distribution function for quarks can be approximated at 
small x using Eqs. (43) and (75), 

xf q/p (x, n) = ^L JdbjdA 0(AV > a 2 ) (87) 

Clearly, the model for H(6, A) contains non-perturbative physics. Furthermore the squared 
wave function 1/A 4 is a perturbative result that should be trusted only for A <C R p . Is 
there any reason to think that Eq. (87) might be reliable at all? 

To examine this issue, first look at the integration range for A. There is a renormalization 
cut A > a/ /I and we may suppose that we consider scale choices such that a/fj, <C R p . The 
integration extends to arbitrarily large A, but once A > 1/Q s (b) we have S(6, A) pa 1 so 
that the integrand is approximately 1/A 4 . This falloff is sufficiently fast that values of A 
greater than 1/Q s (b) are not important in the integration. Now, Q s (b) is proportional to 
the gluon distribution and at small x there are lots of gluons. For this reason, for a central 
impact parameter b, Q s (b) is larger than the normal 300 MeV soft hadronic scale. With 
xG(x) = 10, a s = 0.2 and R p = 4.5 GeV -1 one gets Q s (0, fundamental) pa 0.6 GeV. 3 If we 
were dealing with a large nucleus or with values of x much smaller than 10 -3 , we could have 



3 This value is consistent with the value obtained by comparison with diffractive DIS data in the somewhat 
different approach [8, 10]. 
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quite a lot larger values of xG(x) and thus a larger saturation scale. 4 Additionally, Q s (b) is 
small near the edge of the proton. If we were dealing with a large nucleus, the contribution 
from b near the edge of the nucleus would be less important than for a proton. 

To the extent that Q s (b) is large, the main contributions to xf g / p (x, /i) come from regions 
in the integrations in which the model is anchored in a reliable perturbative expansion. But 
what if Q s (b) is not so large? Then we must face the facts that the model for H(b, A) is 
non-perturbative and that the 1/A 4 squared wave function is a perturbative result applied 
outside the range of validity of the perturbative expansion. We can analyze these problems 
in two ways. First, the dipole interaction with the proton, S(6, A), should be subject to 
scrutiny. Second, we can consider what would happen if we were to work at a higher order of 
perturbation theory. Then we would have new contributions to the partonic state that hits 
the proton, including the possibility that this state contains more than just two partons. 

Recall first the behavior of S(6, A) in the model of Sec. IV C, supposing that the de- 
scription of the incoming partonic state as a dipole with the perturbative 1/A 4 squared 
wave function is exactly right. It is indeed true that S(6, A) cannot be reliably calculated 
perturbatively when A is not small. However, 5(6, A) corresponds to the probability that 
the dipole scatters. When A is large and |6| < Rp, it is likely that the dipole is almost com- 
pletely absorbed, which corresponds to S(6, A) « 1. The model of Sec. IV C for H(6, A) 
has this property. Thus H(6, A) is fixed for small A and for large A as long as 6 is well 
inside the proton. It is certainly true that it is not so well known for intermediate values of 
A and for large or medium A when \b\ ~ R p . In particular, if Q s is not so large one is likely 
to make an error by extending the transparency region to intermediate A. But the effect 
is not dramatic, so that even here there is not too much that one could do to drastically 
change 2(6, A) from the form given by the model. 

Consider now the higher-order states. The original eikonal quark plus an antiquark state 
can become an eikonal quark plus an antiquark plus several gluons, for example. We could 
still define a measure A of the transverse size of this partonic system. The partonic wave 
function would depend on A. It would also depend on other dimensionless shape variables 
that we could call 7. Then we would have a function H 7 (6, A), given by a matrix element of 
multi-eikonal operators, describing the probability for this state, labeled by b, A and internal 
quantum numbers 7, to scatter. This would give an extension of Eq. (87) with the form 

xf q/p (x if j)= fdb fdA 2 ]T^ 7 (A)| 2 ~ 7 (6,A). (88) 
J J 7 

The integral needs renormalization, which can introduce logarithms of A/i. Except for this 
appearance of the renormalization scale /i, the calculation of the wave function |^ 7 (A)| 2 
involves no hadronic distance scales and no masses. For this reason, dimensional analysis 
tells us that |^ 7 (A)| 2 is proportional to A -4 times the logarithms of A/i times dimensionless 
constants and times factors of a s . This suggests, although it certainly does not prove, that 
the squared wave functions |^ 7 (A)| 2 are not larger than the lowest order result. This leaves 
us with the scattering probabilities H 7 (6, A). We do not know the detailed form of these, but 
it is plausible that the complicated states under discussion are almost completely absorbed, 
which corresponds to S 7 (b, A) « 1. This is just the behavior of the simple dipole version of 
the scattering probability, H(6, A), for large A. 



4 Also, Q s (b) is larger if we had a color 8 dipole instead of a color 3 dipole, as would be the case if we were 
to investigate the gluon distribution. See [11] for a recent discussion. 
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These arguments do not establish that Eq. (87) for the quark distribution function must be 
highly accurate if applied to a proton with x ~ 10~ 3 rather than, say, a very large nucleus or 
very much lower x. However, they do suggest that the picture has enough qualitatively right 
features built into it that it should be more useful than would seem from first appearances. 



VI. THE STRUCTURE FUNCTION 

The dipole results for structure functions are known from [4]. The transverse structure 
function Ft is given by [4] 

F ? = ^£^^jf da [l-Ml-*)] JdbjdA 

x | y/ a (l - a) QA K' (y/a(l - a) QA) |* S(6, A), (89) 

where Q 2 is the photon virtuality, and K' Q is the derivative of the modified Bessel function. 
The main difference compared to the case of the quark distribution is that the ultraviolet 
region of small A is now naturally regulated by the physical Q 2 . In Appendix C we sketch a 
derivation of this result along the lines of our derivation for the quark distribution function. 

In the remainder of this section, we relate this formula for Ft to the normal factorized 
form in which Ft is expressed as a sum of perturbatively calculable hard scattering functions 
Ft convoluted with parton distribution functions. For large Q 2 , the integral in Eq. (89) is 
dominated by two integration regions, A ~ 1/Q and A > 1/Q. We discuss each region in 
turn. 

In the case A ~ 1/Q <C R pi one can use the small A perturbative formula for S(b, A), 
Eq. (67). Then the contribution from this region is a certain one loop integral times a s 
times the gluon distribution function. We can recognize that this has the form of a one loop 
contribution to Ft times the gluon distribution. We do not analyze it further. 

The case A ^> 1/Q is more interesting from the point of view of this paper. Let us 
implement the requirement A ^> 1/Q in a crude fashion by inserting a factor 8(QA > c) 
where c is a fixed number of order 1. The only way that we can get a leading contribution to 
the integral for large QA without the Bessel function cutting off the integral is for a(l — a) 
to be small. That is, either a must be small or else 1 — a must be small. We consider the 
case a < 1. To see what this region contributes, we simply neglect a compared to 1 inside 
the integral, 



1 

X A^ 



\V^QA K' (^QA)\' H(6,A). (90) 
Here we can change variables from a to z 2 = aQ 2 A 2 , giving 

1 n r r i r°° 

F T Tq = -E e 'j^ JdbJdA6(QA>c)Z(b,A)— jf dz z\zK' (z)f . (91) 



Using 



dzz\zK' (z)\ 2 = -, (92) 
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this is 



a 

Comparing with Eqs. (43) and (75), we see that we have the lowest order contribution to 
the hard scattering, Ft, times the dipole form of the quark distribution evaluated at a 
renormalization scale of order Q/c. The corresponding 1 - a <C 1 contribution gives the 
same Ft times the antiquark distribution. Thus we see that the leading order factorization 
formula works in the dipole approximation with the quark distribution function defined 
independently according to its definition as the proton matrix element of a certain operator. 5 



VII. CONCLUSIONS 

There is an s-channel approximation for structure functions that is quite standard in the 
literature and is, we believe, well motivated. In this approximation, Ft(x,Q 2 ) is given by 
Eq. (89). This has the form of a dipole scattering probability 5(6, A) convoluted with the 
probability to make the dipole. We have presented a variation of the "saturation" model 
[4, 9] for 5(6, A) in Eqs. (83), (84) and (86). 6 The approximation (89) for F T (x,Q 2 ) seems 
to be quite different from the factorized form applicable at large Q 2 , in which Ft(x, Q 2 ) is 
expressed as a convolution of a hard partonic structure function Ft with parton distribution 
functions. The focus of this paper has been to connect these apparently dissimilar pictures 
by investigating the quark distribution function xf q / p (x,fi) at small x using the s-channel 
picture. 

We have found that the parton distribution function for quarks can be approximated at 
small x using Eqs. (43) and (75), 

xf q/p (x, = ^ Jdb JdA 6(A 2 ^ 2 > a 2 ) (94) 

This has the form of the same dipole scattering function 5(6, A) as in Ft, now convoluted 
with a different probability to make the dipole. In fact, the probability to make the dipole 
is beautifully simple, 

N c 6{A 2 n 2 >a 2 ) 



where a is a calculated number of order 1, Eq. (41), that accomplishes MS renormalization 
for the quark distribution, assuming that [i is large. The power behavior, 1/A 4 , characterizes 
the squared lightcone wave function. 

We have seen not only that the quark distribution has a simple form in this picture, but 
also that the normal lowest order factorized form for Ft relates the dipole expression for Ft 
to the dipole expression for f q / p . Furthermore, the evolution equation for f q / p relates the 
exponent in 5(6, A) to the gluon distribution. 



5 Dipole contributions that arc power suppressed with respect to the leading factorized term are investigated 
in [6] for the Q 2 evolution of the structure function. 

6 The principle refinement is the definition of xG(x,ji), Eq. (72). 
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APPENDIX A: CALCULATION OF u(A) 

In this appendix, we compute the integrals for the function u(A) introduced in Sec. II. 
We begin with Eq. (32), 



( 27 r)Wo 7 P2 7 Pl (A 2 +P1)(A 2 +P?)' ( ] 

We can introduce two Feynman parameter integrals to put the denominators into the expo- 
nent. This enables us to perform the Pj integrals 

u(A) = ^f^l°° dA2 J d2 ' 2e P^J d2 ~ 2e P^ e^~^ Pl -p 2 . 

/"OO /"OO 

x / da.\ exp(-«i(A 2 + p\)) \ da 2 exp(-a 2 (A 2 + pj)) 
Jo Jo 

/•oo rOG roo 

/ dA 2 / dan / da 2 e ~ {ai+a2)A2 
Jo Jo Jo 



(2vr) 6 - 4 ^ 

x 



(~*dA~) / rf2 2tpi ex v(~ ai P 2 + lA ' Pi) 

x ( z ^tJ J d2 ~ 2e P2 exp(-a 2 p 2 - ?A • p 2 ) 

,4e /*oo /* oc /*oo 

/ dA 2 / dan / da 2 e"( ai+Q2 ) A2 
Jo Jo Jo 



4iV c /i 



(27r) 6 - 4e 

9 



x — % 



d 

x I i- 



dA j 
WT 



6-4e 



J d 2 ~ 2t p x exp(-a 1 (p 1 - *A/(2 ai )) 2 - A 2 /(4«0) 
J d 2 -^p 2 exp(-« 2 (p 2 + *A/(2a 2 )) 2 - A 2 /(4a 2 )) 

/*OG /*O0 /*00 

/ dA 2 / / cfa 2 e~ {ai+a2)A2 

Jo Jo Jo 



1-e 



x(-^V-l e exp(-A 2 /(4a 2 )) 



2a 2 / \a 2 J 



/•OO /» 

(4^) 2e A 2 / dA 2 / 
Jo Jo 



' doti 
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1 

— + 






«2 



2 6 vr 4 J J a x J Q a 2 

x exp ( -(ai + a 2 ) A 2 - — [— + — ) ) . (A2) 



At this point, we can perform the A integral 



N c ( A 2\2e a2 [°° dai _ l+e f°° da 2 

u W = 7^=1 w ) A y o — «i / — « 

1 / A 2 / 1 1 



2% 4 v ^ ; i «i io «2 



'2 



x ; exp -— — + — . (A3) 

«1 + «2 V 4 V a l °2, 
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In order to simplify the exponent, we can change variables to $ = l/af. 



P. 



2 Pi 



•2 



P 1 + P 2 ^\ 4 
Now we can change variables to 



J^expf— (& + &)). (A4) 



7 = A + 



r = (A5) 



The inverse transformation is 



The jacobian is 



Pi = n, 

P 2 = (1 - r) 7 . (A6) 



<9(/3i,&) 

= 7- (A7) 



5(7, r) 



Thus 



AT POO pi 

u(A) = ^ (4vr/i 2 ) 2e A 2 ^ 7 d 7 j( dr (ry) 1 "' ((1 - r) 7 ) 



l-e 



2 



1 / A 



2 6 7T 4 

x [ dr (r(l -r)) 1 " 6 . (A£ 



o 



We can perform both integrals with the result 



r(2- e ; 



T(4 - 2e) 

37r 4 A 4 v ; 1 - 2e/3 V ; 



APPENDIX B: AN ALGEBRAIC RELATION FOR EIKONAL OPERATORS 

In this appendix, we seek a relation between the operators Tr(F^F) for the quark and 
the gluon distributions, where F is given in Eq. (45). 

Denote by V and U the eikonal operators in the fundamental and adjoint representation: 

V(z) = F iund .(z) , U(z) = F adj .(z) . (Bl) 
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The following identity holds between V and U at the same point: 

hj a \z) = Tr [t a V{z)t h V\z)] , (B2) 

with t a and t b generators in the fundamental representation. This can be seen by constructing 
the adjoint representation from the product of 3 and 3. 

Using (B2) we can write the trace of two LPs at points x and y as 

Tr [U{x)U\y)\ = U ab (x)U ab (y) 

= 4^.VSK^)tL^(^)^K r (y)4KU^) • (B3) 

Now with the identity 

tijt a pq = ls iq 8 pj - ^-JijSp, (B4) 

we obtain 



Tr[U(x)W{y)] = Tr \V{x)V\y)] Tr [V\x)V(y 

- ( Tr [V{x)V\y)V{y)V\x)] + Tr [V {x)V\x)V {y)V\y)\ } 

+ -^Tr [V(x)V\x)] Tr [V\y)V{y)} 

c 

= Tr [V(x)V^y)] Tr [V*(x)V(y)] - 1 . (B5) 
For the operators that appear in the definition of S, Eq. (44), from Eq. (B5) we get 

^J^Tr [1 - W(A)U(0)] = ^ Re Tr [l - V\A)V(0)} (B6) 



From this general relation we recover the simple ratio (66) in the case of small A, where the 
quadratic term in the right hand side of Eq. (B6) can be neglected. 



APPENDIX C: THE EVOLUTION OPERATOR FOR F T 

In this appendix we derive the s-channel formula (89) for the structure function Ft in 
the same fashion as was done for the quark distribution in Sec. II. 
We start with the definition of Ft, 



Ft — — 

Here 



d 4 ye- l ™(P\Ji(0)Ji(y)\P). (CI) 



P = (P + ,0,0). (C2) 
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Similarly to Sec. II, we rewrite this as 

f d P' + f 

F T = Re(27r)- 3 y — JdP'(P'\U[A}-U[0]\P). (C3) 

Here U[A] is a function of the field operator A, defined by 

U[A] = ^ J dy + J dy x J dy 2 J dy^ J dy^ 6{y^ > yf) 

Xe -*9-v + e -i*P + (y a --«r)(o| Ji( , y", y 2 )J'(y+ yf, yi )|0)^, (C4) 

where the matrix element here is taken in an external potential A. 

Now using the interaction picture with A as the perturbation we have 

U[A] = ^ y dy + y dyi J dy 2 J dy^ J dy 2 Q{y 2 > yf) 

xe-*-»V fa ^foT-vr) ( |t/(oo, y 2 ~) J'(0, y 2 ", y 2 ) 
xU(y 2 ,y^)P(y + ,y^, yi )U(y^,-oo)\0) A . (C5) 

In the approximation that the potential is negligible for large \y~\ while only large positive 
y^ and large negative yf dominate the integrals, this is 



P+ r r r r° f°° 

U ^ ~ 27r J dy+ J dVl J dm J dV ^ J dy ~ 2 



-n y + e -ixP + (y 2 -Vi ) 
x (0| J*(0, yj, y 2 )f/(oo, -oc)^'(y + , yf, yi )|0)^. (C6) 

We understand here that we are going to use the eikonal approximation for U and if we go 
beyond the lowest approximation there will be an effective interval — y^ < y~ < y$ for y~ 
inside the approximation. 

We will evaluate this at the lowest order of perturbation theory for the quantum part 
of the theory. That is, all of the particles are treated as free except for the interaction 
with the external field in A. To carry out this evaluation, we insert intermediate states. 
The intermediate states consist of a quark (momentum k) and an antiquark (momentum p). 
These particles carry spin and color, but we choose a notation that suppresses the spin and 
color indices. Thus we have 

p+ r r r rO roo 

U ^ ^ 2^ J dy+ J dVl J dV2 J dy ^ J dy * 

Xe -iQ-y + e -^P + (v2-yi) (2n)- 12 

x(0|J J (0,y 2 ,y 2 )\p 2 ,p 2 ,k 2 ,k 2 ) 

x (p 2 ,P2, fcj, k 2 \U(oo, -oo)\pi,p u fcf, fei)^ 

x (pi , pi , K , fci I Jj (y + 7 vi , vi) I o) ■ (C7) 
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For the matrix element of U we use the (leading) eikonal approximation, 
(p 2 ,p 2 , k 2 , k 2 \U(oo, -oo)bT ,Pi, ky , fei)^ 

= (27r)>^(pr -P^KKK ~ K) F c (p 1 -p 2 )F(k 1 - k 2 ). (C8) 

Here F is the eikonal factor for the quark, i£ is the eikonal factor for the antiquark. In the 
matrix elements of the current we can use translation invariance to extract the y dependence. 
Then 



U[A] 



P+ f f f f° f°° 
2tT J + J dVl J dV2 J ^ J dy ~ 2 



Here 



xe~ ig y + e - ixP+ (v2 -2/1 ) (27r)~ 10 

x l w-l wf dp >J dh *f dpi J dkl 

x(0\J\0)\p~,p 2 ,k~,k 2 )e- i{p z +k * )y 2 +iiP2+k2) - y2 
xF c (p 1 -p 2 )F(k 1 -k 2 ) 

x(p~,p 1 ,A;-,fc 1 |J i (0)|0)e i(? '" +fc " ) ^ +i( ^ + ^ ) ^- i(pi+fcl) - yi . (C9) 



Pt = A, K = 3z,Pt = K = 3z- (CIO) 



2p-' 1 2k-' r 2 2p~' 2 2k~ 
We can now perform all of the y integrations to get 

x(0| J j (0)\p-,p 2 , k~, k 2 )e~ l ^ +k ^ +liP2+k2) - y2 
xF c {Pi-P2)F{k l -k 2 ) 

x(p-,p 1 ,fc-,fci|J J '(0)|0)e i( ^ +fc " ) ^ + ^ +fe i +)2/ i"- i(pi+fel) - !/1 



xP + +p 2 +k£ xP + +pf + kf 
x5(p 2 + k 2 ) 6{ Pl + fci) 8(q- - p~ - k~). (Cll) 

For the minus-momenta we write p~ = aq~ and k~ — (1 — a)q~ . With the use of the delta 
functions, the momenta are 

( P 2 

h = Ui-V ,(1 ' a)9 "'" pi 

( Pl - 
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This gives 

x g— ^ i— r+ p+ ~' + - , + Tr{F c ( Pl - p 2 )F(p 2 - Pl )} 

a;P+ + p 2 + k 2 xP + + + k{ 

x(0\P(0)\p 2 ,k 2 }(p 1 ,k 1 \P(0)\0}. (C13) 
For the matrix elements of J J , we can write 
(0\P(0)\P2,k 2 )(puh\P(0)\0) 

= 5Z e « 5 s lS2 5 s ' lS 'u(k2,S2)i 3 v(p2,s' 2 ) v(pi,s\)^ J u(k u si). (C14) 

Now we can insert 

$ s ' lS ' 2 = v(p 2 ,8 , 2 )'y-v(p 1 ,s' 1 )/(2p-). (C15) 

This leads to 

(0|^'(0)|p 2 ,fc 2 )( Pl , fci|J'(0)|0) 
= 4a(l -a)( g -y E e >{^ 2 7^i7 J >i7-fe} 

= 4 E e ' 1 ~n a(1 T ) pi ' P2 - ( C16 ) 

all — a) 

Thus 

rrr ,, 1 2 2P+ y 1 , /\ /\ l-2a(l-a) 

U W - ^ta^^l d* J dp 2 Jd Pl a2(i _ a)2 Pl - P2 

—i —i 
xP + +p% + k£ xP+ +pf + kf 
xTr{F c ( Pl - P2 )F( P2 - Pl )}. (C17) 

We can rewrite the energy denominators to obtain 

U W ~ ~^E e * da J dp2 S dpi [ 1 - 2a (i-«)]pi-P2 

1 1 

X a(l -a)Q 2 + P 2 2 a{l-a)Q 2 + P \ 

xTr{F c ( Pl -p 2 )F(p 2 - Pl )}. (C18) 

With 

F c ( Pl - p 2 )F( P2 - Pl ) = J dbdA e i(Pl-P2)-(6+A/2) e i( P2 -pi)-(6-A/2) 

xF c (6 + A/2)F(b- A/2), (C19) 
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we get 



xTr{F c (6 + A/2) F(6 - A/2)} 



7 «(i 



P2 



x / e^ A — P l 2 . (C20) 

We now take the hadron matrix element (C3), and use the definition (76). Then performing 
the integrations over pi and p 2 gives Eq. (89). 
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